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Abstract 

One approximation is made to describe a M+l electron many-body wavefunc- 
tion by a M electron many-body wavefunction and a single electron wavefunc- 
tion. Under this approximation, we have derived the Coulomb energy which 
relates the exciton energy E^xc i^^ quantum dot with the quasiparticle band 
gap (defined as the difference between the ionization energy and the electron 
affinity), and the Coulomb energy which relates E^^xc with the single particle 
eigen values. We found that these two Coulomb energies are different. We 
have compared our results with the formulae used in different groups, which 
are proposed either ad hov, or from the classical electrostatic point of view. 
We found important difference between our results and the classical formu- 
lae. Finally, under the above approximation, we provided an effective single 
particle Hamiltonian, which gives the quasiparticle band structure in a bulk 
system. 
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I. Introduction 

The screened Coulomb interaction in a bulk semiconductor between an excited electron 
and a hole has been studied theoretically 30 years ago [|T],3 from many body point of views. 
As a result of such theoretical analysis, the screening of the Coulomb interaction can be 
expressed using the dielectric response function of the bulk system. Recently, nanometer 
scale finite system has attracted a lot of research interests . The confined excited electrons 
and holes within the finite systems have enhanced Coulomb interactions comparing to the 
bulk material [Q. However, there are different opinions about how to screen the Coulomb 
interactions in various cases. While the screening of the Coulomb interaction is clear when 
the electron is outside the nanostructure, it becomes clouded from a classical point of view 
when it enters the nanostructure. One problem is that the classical dielectric screening effects 
can already been partially represented by the single particle eigen energies. Thus combining 
the classical electrostatic model with the single particle Hamiltonian does not always provide 
a clear picture free from ambiguities. The more rigorous approach is to derive the single 
particle effective Hamiltonian and the screened Coulomb interaction from the many body 
Hamiltonian. However, the Feyman diagram technique used in the bulk study thirty years 
ago can not be directly applied to the finite system without complicated modifications of the 
Green's functions according to the new quantum dot boundary conditions. In this article, we 
will present a derivation of the screened electron-hole Coulomb interaction in a nanosystem 
started from a many body Hamiltonian. Instead of basing on the Green's functions, our 
derivation is based on many particle wavef unctions. We will restrict ourselves to the cases 
that the excited electrons and holes are strongly confined by a dimension nanostructure 
(quantum dot). That means, the correlation between the excited electron and the hole can 
be ignored |^,^ 

One of the earlier works on the quantum dot Coulomb interaction and its screening 
effects is given by L.E. Brus ||^. Using a classical model of the dielectric screening and a 
single particle Hamiltonian, Brus derived the change of electron affinity of a quantum dot 
(spherical, with radius R) relative to the bulk system: 
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En - En+1 = E^/f - eg - ^-P'- (1) 

Here Em is the total energy of the quantum dot with M electrons. E^ is the neutral 
quantum dot. Thus Ej^ — -Eat+i is the electron affinity of the quantum dot. E^f^ is the 
electron affinity of the bulk, eg is a single electron confinement energy in conduction band, 
which equals 2m*B? u^'i^^ effective mass model of electron effective mass m*. i^P' is an 
electrostatic energy representing the interaction between the electron and its image charge 
produced by the quantum dot medium with dielectric constant e. Similarly, for ionization 
energy, he has 



En-i -En = E^i!^ + e, + ^P'. (2) 



Here eh is the single hole confinement energy, which equals 2m* under a single band effective 

h 

mass model of hole effective mass m^. The last term in Eq(2) is the same as the last term 
in Eq(l). Then, using the conventional definition, the quasiparticle band gap E'^^ equals the 
difference between the ionization energy the the electron affinity: 



Ef = £^^+1 + En-1 - 2En = E';''' + ee + e, + ^P' . (3) 

Here, we have used the fact that the bulk band gap E^^^^^ equals the difference between bulk 
ionization energy E\'^ and the bulk electron affinity -^a/'/- In this paper, we discuss only 
the cases where the quantum dot is surrounded by vacuum. Then if the dielectric constant 
of the quantum dot e is much larger than 1, we have (following Ref. PJ): 



V = (1-1)1 + ^(1^). (4) 

In a second paper 0, Brus presented the result for the exciton energy E^xc in a quantum 
dot (which equals the photon energy needed to excite an electron from the valence band to 
the conduction band). 
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where — represents the direct Coulomb integral between the electron charge and hole 
charge, screened by e, using the effective mass wavefunctions. The last term in Eq(5) is said 
to be the time averaged instantaneous dielectric energy while the classical particle moves 
around inside the quantum dot. Comparing Eq(3) to Eq(5), with the help of Eq(4), we have 

E,,, = Ef - — -il--)-. (6) 

Note that, classically, the last term in Eq(6) represents the Coulomb interaction energy 
between a spherical charge density inside the quantum dot with its induced surface charge. 
The total of the last two terms in Eq(6) represents the electrostatic interaction between 
the confined bare electron charge density with a screened hole charge (i.e, the total charge 
density of the -Eiv-i system). 

However, the above classical picture is not universally accepted. For example, in a recent 



paper of Ogut, etal [ly], the connection between i?^^ and E^xc is not given by Eq(6), instead 
a simple Coulomb interaction — is used (in the effective mass limit): 

E.XC = - 1|. (7) 

The difference between Eq(7) and Eq(6) is very large (up to 2 eV for R ~ 7A), because 
there is an unscreened Coulomb interaction —1/R in the last term of Eq(6). Using Eq(7), 
Ogut etal obtained an exciton energy Eg^c much larger than the one obtained from the single 
particle empirical pseudopotential method (EPM) |TlI] (which is yield via Eq(5) but without 
its last term, and with eg, calculated via single particle EPM Hamiltonian) . Based on this 
difference, they concluded that the empirical pseudopotential calculation is wrong because 
a size dependence of the GW selfenergy J2 is ignored. However, if Eq(6) instead of Eq(7) is 
used, one finds that the result from Ref. [|10| is almost the same as the EPM result. Then 
their conclusion will be false and EPM calculation will be all right. 



To settle the dispute between Eq(7) and Eq(6), and to test the other classical formulae 
[e.g, Eq(5)], it will be very useful to derive the above relationships from many body Hamil- 
tonians. This will be provided in the rest of this paper. We found out that Eq(6) is correct 
instead of Eq(7), and Eq(5) is correct after deleting its last term. 



In this section, we will derive our basic formula to be used in later sections. This formula 
is based on an approximation that the many-body wavefunction $m+i of a M+l electron 
system can be separated into one single particle wavefunction [i.e, ip^ (top of valence band 
state) or ipc (bottom of conduction band state)] and the rest of the M particle wavefunction 
$M- More specifically, we have: 

^M+l{xi, X2, Xm, Xm+i) = \^Mlp >= , [^Mixi, X2, Xm)'4'{xm+i) 



Here, we have used |$m'0 > to denote the break down of the total wavefunction $m+i 
into $M and ip. x = {r,a), and r is the three dimensional Cartesian coordinate, and a 
is the spin index. Note that through the definition of Eq(8), ^m+i{xi,X2, ...,xm,xm+i) 
is antisymmetric, provided $m is antisymmetric. In general, ^m{xi,X2, ■■■,xm) can be 
represented by a linear combination of different electron configurations (i.e, single slater 
determinates from a single particle orthogonal basis set including ip) . Any slater determinate 
which consists of ip will be eliminated in Eq(8) by the antisymmetry operation. Thus, it 
does not lose any generacity (in terms of variational degree of freedom for $m+i) to exclude 
the single particle orbital ip from the configurations of ^m{xi,X2, ...,xm)- In another word, 
we can restrict ^Mixi,X2, ...,xm) to be orthogonal to ip{x), i.e. 



II. The basic Formalism 



^m{xm+i, X2, ■-, xm)iP{xi)... - ^m{xi, X2, xm+i)iP{xm)]- 



(8) 




(9) 
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This orthogonal condition is also called strong orthogonal condition [^,0. Notice that, 
since we have not restrict any degree of freedom in $Af+i of Eq(8) by introducing Eq(9), 
Eq(9) should not be considered as an additional approximation. The approximation, if 
any, has already been made in Eq(8). By using Eq(8), we have neglected the correlation 
between the single particle state ip and the rest of the system This is one approximate 
way to define an quasiparticle wavefunction ip{xi) out of a many particle system. The 
approximation of Eq(8) has been called group function approximation [1T2| , |T3| , and has been 
used to study the effects (e.g, screening) of one group to another group. Thus, it is natural 
here to use it to describe the Coulomb screening effects of the background system ($a/) 
to the quasiparticle (and exciton) system {jp). Note that, Eq(8) goes beyond Hartree-Fock 
1^ approximation by retaining correlations within $jv/- Since our focus is to study the 



dielectric screening of the Coulomb interaction, which can already be described by Hartree 
Fock approximation, we expect Eq(8) to be adequate to represent the main physics of 
interest here. Similar approximations of Eq(8) has been used in the derivation of Extended 
Koopmans' Theorem by Day, Smith, and Garrod []T3|, by Morrell, Parr and Levy |]TB[, and 
by Kent, etal [|r^. In Eq(8), one gets the wavefunction $a/+i of the M+1 electron system 
by adding one electron ip to the M electron wavefunction $Af- Similary, one can get the M 
electron system wavefunction $a/ by elminating ip{xi) from the M+1 electron wavefunction 



<l>Af (Xi, X2, XAf) = j ij*{xM+l)^M+l{Xl, ^2, Xa/, XM+l)d^XM+l- (10) 

Using the orthogonal condition Eq(9), one can derive Eq(lO) from Eq(8). Equation (10) is 
also used in the derivation of the Extended Koopmans' Theorem ||l5Hl7|. 

To appreciate the approximation embodied in Eq(8), we can take a look at the density 
matrix of the M+1 system. Under Eqs(8) and (9), Pm+i{.x,x') = ip{x)%l)*{x') + pMix,x'), 
here pm+i and pa/ are the density matrix of the M+1 and M electron systems respectively. 
Besides, following the orthogonal condition of Eq(9), / pm{x,x')iP{x')cPx' = 0. As a result. 
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ipix) is a natural orbital of pm+i with an occupation number of 1. As well known [|T^], the 
occupation numbers of the natural orbitals are usually less than 1, thus the Eq(8) is definitely 
an approximation. However, this approximation is similar to the approximations which use 
one single particle wavefunction to represent a quasiparticle in a many-body system. These 
include: the Hartree-Fock theory [|T1|, empirical pseudopotential theory |TT|, the Kohn-Sham 



representation of the density functional theory |20| and in some extent the GW theory JTS 



To include the correlation between ip and $Af in Eq(8), one has to add in Eq(8) other 
configurations I^'mV'' >• However, if the quasiparticle is represented by a unique single 
particle wavefunction ip, then ip' must be ip. As a result, |$Af^ > +\^'m'^' >= I'^X/^ >5 
where = + ^'m^ Eq(8) remains unchanged. 

The Hamiltonian for the M+1 particle system can be written as 

M+1 M+1 M+1 2 

Hm+1 = E V.^ + E yUn) + E TT—r-r W 
^ i=i i=i ijtj l'« ' i\ 

Using Eq(ll) and Eq(8), we can evaluate < ^m+i\Hm+i\^m+i >■ After some algebras, it 
can be expressed as 

< ^m+i\Hm+i\^m+i >= E^>+ < ^m\Hm\^m > 

+ I PM{r)Vt^i{r)d'r - j Kt(^, x')-pM{x, x')d'xd'x', (12) 
where Hm is just the Hm+i of Eq(ll), but change the M + 1 to M. i?^ is defined as 



2 

V^^i{r) is the Coulomb potential due to ip 



i)*{x)d^x. (13) 



VLi{r)= I -^Mx')\'d^x\ (14) 



PM^r) is the total charge density of $ 



M, 
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pM{r) = \^M{{r,cr),X2,...XM)\^d^X2...d^XM, (15) 

V^{x,x') is the exchange nonlocal potential due to 



Kt{x,x') = y-^^r{xmx'), (16) 

and pm{x,x') is the density matrix of $m, 

Pm{,X,x) = M j (^Mix,X2, .■,Xm)^*m{x',X2, ■.,XM)d^X2...d^XM- (17) 

Note that pM{r) = Y.aPM{ra,ra). 



Equation (12) has been derived before as in the group function theory Under the 
assumption of the wavefunction ^^z+i of Eq(8), one can variationally minimize the total 
energy of Eq(12) with regard to $Af and if) under the constraint of Eq(9). Doing so, both 
$A/ and ip can be solved selfconsistently from Eq(12). However, symmetries might need to 
be used to keep the ip to be the desired single particle wavefunctions (i.e, the top of valence 
state and bottom of conduction state). Another approach is to calculate ip from single 
particle Hamiltonians, e.g, quasi-particle GW method [|18[, or empirical pseudopotential 
method ||lT|. Later in Section IV, we will provide a new single particle Hamiltonian Hs 
[Eq(40)] to calculate ip. For a given ip, we can minimize Eq(12) with regard to "I'm, thus get 
a many particle equation for $m under the perturbation of Coulomb potential Vf^^i{r) and 
the exchange potential V^{x,x'). 



III. Coulomb energy associated with the ionization energy and electron affinity 

The exciton energy E^^c can be directly defined as E^f — E^^. Here, E^ is the N electron 
ground state energy and E"^ is the total energy of the N electron system which contains 
one exciton (i.e, an conduction band single particle state ipc is occupied and an valence 
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band single particle state ipv is unoccupied). Then the Coulomb energy which connects E^p 
(= En+1 + En-i - 2Em) with Ee.xc is 

Ecoui = Eexc — E'^P = [El, — Em^i) — {Enj^i — En), (18) 

where, E^-i and E^^i are the ground state energies of N-l and N+1 electron systems. 
Note that in E^^, E^, Ejq-i and Ejq+i systems, the ionic potential Vion^r) in Eq(ll) is 
the same. The difference is the electron occupations. Using the approach outlined in the 
previous section, we will approximate the wavefunctions of E*j^, Ejy, E^-i and E^+i systems 
as \^%_iipc >, I'^'at >, I'^'Af-i > and l^l^ipc > respectively. Using the orthogonal conditions 
of Eq(9), and applying it equally to all the systems in Eq(18), we have and 
$^ all orthogonal to ipc [Eq(9), replacing ip by ipc, and $a/ by these four wavefunctions]. 
In addition, and should be orthogonal to ip^ [Eq(9), replacing ip by ipv], so 

that a hole exists in these systems. Notice that, we have used the superscription "0" to 
indicate that the corresponding wavefunction is a minimum energy variational solution 
of < >. Similarly, the superscription "c" indicates that the corresponding 

wavefunction is the variational solution of Eq(12) under the perturbation of tpc [through 
V^^^^i and V^"]. Notice that, only in the strong confinement quantum dot size region 0, we 
can write $^ as l^'j^^^ipc > through Eq(8). In this size region, the electron and the hole 
are not correlated. As a result, we don't have to write $^ as a summation of configurations 
l^N-i'^c > with different electron [tpc) and hole (tp^j) wavefunctions. After the definitions of 
the wavefunctions, we can replace E^^, E^, E^-i and -Eiv+i by < ^^^^^-j/'d-f^Af I'^'at-iV'c >, 
< <^%\Hn\^1] >, < ^Vil^^-il^^w-i > and < $^?/^c|i^7v+i|$^'0c >, respectively. Using 
Eq(12), we have: 

El - En-i = E^^+ < ^%_^\HN-i\^%^i > - < <!>%_i\Hn-i\^%_i > 

+ J PN-i{r)Vti{r)dh- - J p%_,{x, x')Kf{x, x')d'xd^x' , (19) 
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and 



Em,. -En = E., + < ^Ij\Hn\'^m > - < ^^m\Hn\^Ij > 



+ 1 PN{r)Vti{r)dh - I p%{x, x')V,f{x, x')d'xd'x'. (20) 

The changes from ^% to can be described in first order by perturbation 

theory, under the external perturbation potentials V^^i{r) and V^" in Eq(12). For a N-1 
electron system, the energy changes from < ^^j^_^\Hn_i\^%_^ > to < JifAr_i|$^_^ > 
in response to these perturbative potentials. The amplitude of this energy change is in 
the same order of the Coulomb interaction of the perturbation charge (V'c) with the re- 
sponse charge of the system (which is also in the order of ipc)- other words, AEn_i =< 
$5v-il-f^iv-i|^Ar-i > — < ^n-i\Hn-i\^n-i > scales as 1, instead of N. This is also evident 
from Eq(19). Note that E^ — E^-i scales as 1. At the right hand side of Eq(19), combining 
the E^^ term defined in Eq(13) with the / p'N-i{r)V^J^i{r)d^r- J p%_^{x, x')V^j'{x, x')(fxd^x' 
term, one gets an energy similar to that of the Hartree Fock single particle eigen energy. 
Thus, this combined term also scales as 1. As a result, the AE'^v-i scales as 1. The 
same is true for AE'at =< ^%\Hn\(^% > - < <l>^|i/Ar|<l>^ >. The difference between 
AEn_i and AE^ is the difference between N-1 and N electron system, which should scale 
as AEn/N oc 1/N. Thus, this difference is negligible comparing to the Coulomb energy 
Ecoui which we are calculating. As a result, subtract Eq(20) from Eq(19), we have 

Ecoui^ j[p'N-i{r)-p'N{r)]Vti{r)d'T 

- j [p^_i(x, x') - p%{x, x')] V,f{x, x')d'xd'x'. (21) 

Further more, following the same perturbation argument, we have p$f{r) —p%{r) ~ Pat-i (^) ^ 
PAr-i(^) ^-iid pj^{x,x') — p%{x,x') ~ p%_i{x,x') — p%_i{x,x'). Again, the approximation has 
an relative error oc thus negligible. Then, we can change Eq(21) to 

Ecoui^ J[p%-i{r)-p%{r)]Vt,{r)d\-J [ff^_,{x,x')-ff^{x,x')]v,f{x,x')d'xd'x'. (22) 
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This is our central result. Notice that, p%_i{r) is the total electron charge density of a N-1 
electron system (containing one hole), and p%{r) is the total electron charge density of a N 
electron neutral system. Thus, p'^J^{r) = p%_i{r) — p%{r) is the screened hole charge (the 
bare hole charge plus the response screening charge) in the quantum dot. Note that 



pl''{r)d'r = -1. (23) 

Thus, if we ignore the exchange interaction in Eq(22) (which is usually much smaller than 
the Hartree interaction), and express V^l^i^r) using Eq(14), we have 

i?e.c = + / ^^^^^^^i^ciV.ciV, (24) 

J \ri — 

where, pdr) = J2a \4'c{i^'^)\'^- Iii the case where phenomenological dielectric constant e can 
be used to describe pl^'^{ri), we found that Eq(24) is in agreement with Eq(6), not Eq(7). 
More, explicitly, in that case we have 

p^'^(r) = -^p^(r) for r < R 

= -(1 - -)6{r - R) /A-nR^ for r R 
e 

= for r> R, (25) 

where R is the quantum dot radius and Pv{r) = J2cr\i^v{ra)\'^. Using the effective mass 
charge density Pv{r) = Pc('") = sm^(7rr/i?)/27ri?r^, we yield Eq(6) from Eq(24). 

In the empirical electrostatic derivation of Eqs(l)-(2), it is not so clear whether the 
empirical single particle eigenvalues eg, eh [in Eqs(l)-(2)] already include the effects of elec- 
trostatic energies. For example, when one electron is removed from the quantum dot, how 
to relate the single electron eigenvalue and classical electrostatic energy to the change of the 
total energy in the system is problematic 0J^. Here, deriving Eq(24) from the many-body 
wavef unctions, we do not have all these conceptional difficulties. That is the biggest advan- 
tage of the current derivation comparing to the classical empirical derivations 0J^. Notice 
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that, the dielectric constant e is not derived here, unhke in the case of bulk exciton screening 
0. This, of course, doesn't mean that we have no screening in our formalism. The dielectric 
screening of the Coulomb interaction is enclosed in pl^'^ of Eq(24) in our formula. We just 
didn't derive the detail screening function in our formula (i.e, how to calculate p^^*^ from p^ 
by a detail e model). That is not the focus of our current study. In the cases of Eqs(6) and 
(7), the simple response model of Eq(25) is assumed. It is in this context, e is used here, 
and Eq(6) is found to be correct. 



IV. Coulomb energy associated with the single particle eigenvalues 

Had proved Eq(6), now we like to test Eq(5). What important here is to define an 
effective single particle Hamiltonian, for which eg and are its conduction band minimum 
and valence band maximum eigen energies. Let's start from the definition of the exciton 
energy: E^xc = E*^ — E^. Now, we will rewrite the N particle ground state wavefunction 
<l>^ of En as \^n-i'^v >, with $5^- 

being orthogonal to both ipy and ipc- Using Eq(12), En 

can be written as: 

En — E^^ =< ^^j^_^\Hn-i\^n-i > 
+ / PN-iir)Ktiir)d'r - J ff^_,{x,x')Vfj{x,x')d'xd'x' . (26) 

Here, E^^ is evaluated from Eq(13). Similarly, like before, the exciton wavefunction can 
be expressed as \^'j^_iipc >, with being orthogonal to ipy and ipc- Then, E^^ can be 

written as: 

E% - E^^ =< ^% _;^^\Hn -11^% > 

+ J P%-i{r)Kti{r)d\ - J p%^,{x,x')Vf,^{x,x')d'xd'x'. (27) 

Comparing Eq(27) with Eq(26), we find that we can obtain ^n-i from ^^n-i by applying 
potentials V^tli - Vlli and Vf^^ - Vf^^ into Eq(26). More specifically we can define 
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E{(3) =< $^_Jiff_,^|<_, >= 

< K-,\Hn-i + [(1 - P)V^:i + (3Vt^ - [(1 - + /3KtlK-i > ■ (28) 

In Eq(28), is defined as the minimum energy variational solution of the Hamiltonian 

H^ff, while subjected to orthogonal condition Eq(9) to V'c and '^y. Using Eq(26) and Eq(27), 
we have E{I3 = 0) = En — E^^ and E{(5 = 1) = E*j^ — E^^. Now, using the "adiabatic 
integration technique" , we have 



E{P = 1) = E{f3 = 0) + /' 

Jo 



dE{(3) 
d(5 



dp. 



Note, in Eq(28), ^^r-i is the variational solution of -f^f//, as a result. 



(29) 



< 



N-l 



dp 



|-"e//I^JV-l >-< y^JV-l|-"e//l gp 



>= 0. 



Then, we have 



E{P=1) = E{P = 0)+ <<_^|^^|<_^>rf/3 

Jo 



dp 



yipv 
coul 



(30) 



- j -p%_,{x, x') [Ktl^, x') - J V,f{x, x')] d'xd'x'}dp. (31) 

Pn-1 and p^-i are the density and density matrix of the N-l electron system under the 
perturbation of P[V^^li — V^t«/] and PlVf^" — V^"] (the nonperturbated values are and 
Piv-i)- These perturbative potentials are due to charge density /3[|'0c(3:^) P — |V'i)(^)P] and 
density matrix P[il;c{x)tlj*{x') — tljv{x)'ilj*{x')]. Here, we will only consider cases where a 
phenomenological macroscopic dielectric constant e can be used to describe the response 
of the N-l electron system to these perturbative charges. These perturbative charges can 
be considered as external to the N-l electron system. Thus, under this phenomenological 
description, the N-l electron charges can be expressed as 
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P^-i(r) = pl-,{r)+P{- - l)E[l^c(r,a)r - |^.(r,a)p] (32) 
-p%_,{x,x') = pl_^{x,x')+P{-^ - l)[Mx)re{x') - M^mix')] (33) 

Note that, p^^iir) — J2a Pn-i[{''^t'^)j ('^j'^)]) thus the dielectric constant e used to describe 
the response of p^_i(r) [Eq(32)] and p%_i{x,x') [Eq(33)] must be the same. Also notice 
that, the dielectric constant e used in Eqs(32),(33) is consistent with the conventional def- 
inition of e, and there is no surface response charge as in Eq(25) because the net charge of 
X!^[|'0c(?', cr)P — |'0i;(?', cr)P] is zero. Now, substitute Eqs(32)-(33) into Eq(31), carry out the 
/? integration, we have: 

E*M — = Em — E^^ 
+ / {pl-.{r) + \{\-l)Y.\mr.^)?-\Ur.^)\%Vt^^^^ 

+ / {f^_,{x,x')-r\{\-mc{x)rc{x'y^^ (34) 

Note that p%,{r) = + \i>v{r,a)\'^ and p%{x,x') = p"^_^{x,x') +il!y{x)'il!*{x'), where 

p%{r) and p%{x,x') are the neutral system N electron ground state charge density and 
density matrix. Using this relations, and E^xc — E% — E^, we can derive from Eq(34) that: 



Ee.c = E^-E^^ + j pl{r)[Vti{r)-V^^i{r)]d\- j pl{x,x')[Vt{x,x')-V^j{x,x')](fxSx' 

_ 1 1 \Ux)?mx')\\.^^.^, + 1 /• ^c(x)v..(x)v.:(x>:(xo^3,^3,. (35) 

e J \r — r'\ e J \r — r'\ 

Similar to Eqs(14) and (16), we can now define Coulomb and exchange potentials due to 
the N electron ground state charge density p%{r) and density matrix p%{x,x'): 

Vliir) = J y^^P%ir')d'r' (36) 
K''Ax,x') = y^^f^M{x,x') (37) 
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Then, we can defined a effective single particle Hamiltonian Hg'. 

=< ^{Hsl^p >= 

J r{x){6{x - x')[-^V' + VUr) + V^Lir)] - VfAx,x')}^{x')d'xd'x'. (38) 
Using Eq(38), and Eqs(14),(16), the first part of Eq(35) can be simplified, and it leads to 

e J \r — r \ 

+1 / ^c(x)^.(x)^:(xo^:(xo^3^^3^. (39) 

e J \r — r'\ 

where e^^ and are the eigenvalues of the bottom of conduction band and top of valence 
band of the effective single particle Hamiltonian Hg in Eq(38). This concludes our major 
result for this section. 

Comparing Eq(39) to Eq(5), we notice that: (1) The bulk band gap E'^"''^ in Eq(5) has 
been absorbed in the definition of the single particle eigen values e^^ and e^^. (2) The last 
term in Eq(5), which represents the classical instantaneous dielectric energy, does not exist 
in the current result. It is thus quite plausible that this term should not exist in Eq(4) either, 
which affects the electron affinity [Eq(l)], ionization energy [Eq(2)] and quasi-particle energy 
[Eq(3)]. 

Equation (39) confirms the conventional way to calculate the exciton energy via the 
single particle eigen values, e.g, as in the empirical pseudopotential approach ||lT|. Both the 



Hartree Coulomb interaction and the exchange interaction exist in Eq(39). Interestingly, 
following the assumption of Eqs(32)-(33), the exchange interaction is screened in Eq(39) just 
like the Coulomb interaction. This might shed some light on the long standing controversy 
about whether the exchange interaction should be screened [p],p|,p!9| [notice also the second 
term in Eq(22)]. Also important here is that we have provided a definition of an effective 
single particle Hamiltonian Hg in Eq(38) using the N electron charge density and density 
matrix. Equation (38) is like the variational equation obtained from Eq(12) for a given 
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$Af- The difference is that, here, the same $Ar is used for both ifj^ and ijjc- This Hs 



can be compared with the conventional single particle Hamiltonians e.g, EPM |]TT|, local 
density approximation (LDA) |^, quasi-particle GW calculation [|I^ and Hartree-Fock 
equation |T^. Hg is almost the same as the Hartree-Fock equation, but that the uncorrelated 
HF exchange potential V^fjp{x,x') = Y.i(^occi'i{^)'^i {^') ^as been changed to the 
correlated exchange potential Vj^{x,x') of Eq(37). According to Eq(39), the eigenvalues of 
Hs should provide the band structure (at least the band gap) of a bulk system. This remains 
to be tested. 

Equations (6), (24) and Eq(39) represents two different approachs to calculate the exciton 
energy in a nanostructure. In the paper of Ogut etal the first approach is taken. 
Unfortunately, Eq(7) instead of Eq(6) was used. If Eq(6) were used, the result of Ogut etal 
would be very close (within ~ 0.1 eV) to the results obtained via EPM calculations |jll[ and 
Eq(39) (a constant 0.68 eV Si LDA band gap correction needs to be added to E'^J'j^j^^^R) in 
order to get this good agreement). 



V. Numerical test of the formulae 



Using local density approximation method, Ogut, etal ||T0[ have calculated the E^+i + 
i?iv_i — 2Em and the single particle Kohn-Sham eigenvalues. If we approximate the Kohn- 
Sham eigenvalue as the single particle energies e^X= ^c) and e^„(= e„) given by Eq(38), 
then we have an good case to test our equations (24) and (39). According to Eqs(24) and 
(39), we have (ignoring the exchange interactions) 



S = {En+i + En^i - 2En) - (ec - e.) = (1 - -T5^)n (^0) 

t\ri ) ri 



We have used the same symbol S in Eq(40) as in Ref. ||T0|, although we do not mean that 



it is the GW "selfenergy" as claimed in Ref. ||T0|. According to Eq(40), S is simply an 



electrostatic energy between the electron charge and its induced surface charge. Although 
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there is a well know LDA error for the band gap energy, this error exists in both {En+i + 
Eiy_i — 2En) and {e^ — e^), thus should be cancelled. As a result, Eq(40) should still be 
valid for LDA calculations. 

In Fig.l, the quantity {E^+i + E^-i — 2i?Af) — (ec — e^) is plotted as a function of 1/R. 
This quantity is compared with (1 — -^^)-^- The agreement is quite good. Here, e{R) is a 
function of the quantum dot radius R. This function is calculated in Ref. [^], and can be 
expressed as e(-R) = 1 + (11.4 — 1)/(1 + (a/i?)'). Here, we have used e(i?) which corresponds 
to the total polarizibility of the quantum dot (e^ in Ref. l2ll]), i.e, a = 4.25A and / = 1.25. Of 
course, the agreement of Eq(40) depends on how good is the phenomenological macroscopic 
description of the quantum dot screening [Eq(25)]. There is no reason to believe that Eq(25) 
is exact for a small quantum dot. After considered all these uncertainties, the agreement in 
Fig.l is quite good. 

Notice that, the LDA Kohn-Sham single particle band gap ec — is almost the same 
(within 0.1 eV) as the single particle band gap calculated from empirical pseudopotential 
after a 0.68 eV band gap correction is added to LDA result. This can be confirmed by taken 



the data from Ref. |TD| and Ref. |TT| . It has also been confirmed separately by Delley etal in 
Ref. [^. Then, using the good agreement between LDA {Ejy+i + Ejy^i — 2Ejq) — (cc — e^) 
and (1 — l/e(-R))/-R as shown in Fig.l, we know that the exciton energy E^xc calculated from 
LDA quasiparticle energy through Eq(6) [plus 0.68 eV correction], should be the same as 
Ef^xc calculated from empirical pseudopotential single particle eigen values through Eq(40) 
ilTI. 



VI. Conclusions 

One approximation [Eq(8)] to the many-body wavefunction is presented. This approx- 
imation allows us to define a single particle wavefunction in a many-body system. This 
approximation is more accurate than the Hartree-Fock Slater determinate for the many- 
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body wavefunction. Under this approximation, we have derived the Coulomb energy needed 
to relate the quasiparticle energy E^^^ to the exciton energy Ef.xc (= the optical transition 
energy) [Eqs(22),(24)]. In the limit where the phenomenological description of the dielectric 
screening is valid, we found that the correct formula for this Coulomb energy is Eq(6), not 
Eq(7). We also derived the Coulomb energy which relates the single particle eigen values 
with the Eexc [Eq(39)]. We found that the classical instantaneous dielectric energy in the 
last term of Eq(5) does not exist in our currently derived formula [Eq(39)]. Under the 
assumption of Eqs(32)-(33), we found that the exchange interaction in Eq(39) is screened 
as for the Hartree interaction. Using Eq(24) and Eq(39) respectively, we found that the 
optical transition energy obtained from the LDA E'^^ (plus the LDA band gap correction), 
is almost the same as the result obtained from the EPM single particle eigen values. Thus, 
the conclusion made in Ref. [|1^], that EPM misses the change of selfenergy with size R, is 
incorrect. Finally, we presented an effective single particle Hamiltonian Hs [Eq(38)], which 
under the assumption of Eq(8), provides the band structure of a bulk system. Further 
testing of Eq(38) is needed. 

I thank my colleagues A. Zunger, A. William, A. Franceschetti and R. Needs for useful 
discussions. This work was supported by the Office of Energy Research, Material Science 
Division, U.S. Department of Energy, under Grant No. DE-AC02-83CH10093. 
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FIGURES 

FIG. 1. The LDA calculated (-Eat+i + ^at-i — 2En) — (cc — e^) compared with the surface 
polarization Coulomb interaction (1 — l/e{R))/R. 
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